Disorder-induced critical phenomena 
— new universality classes in Anderson localization — 
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The Anderson metal-insulator transition is a continuous phase transition driven by disorder. 
It remains a challenging problem to theoretically determine universal critical properties at the 
transition. The Anderson transition in a model with a discrete sublattice or particle-hole symmetry 
belongs to one of seven universality classes which are different from the three well-known standard 
ones. Here we review our recent theoretical work on these new universality classes in (quasi) one 
and two dimensions. 
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I. ANDERSON LOCALIZATION AND THE 
DYSON SINGULARITY 

According to Landau's Fermi-liquid theory, low-energy 
electronic excitations in a solid are weakly interacting 
quasiparticlcs that are adiabatically connected to excita- 
tions in a free electron gas. By combining Bloch's the- 
orem with the Fermi statistics obeyed by electrons, one 
knows if the crystalline state is a perfect metal (having 
an infinite conductivity) or a band-gap insulator at tem- 
peratures low enough that the electron gas is degenerate. 

It is an empirical fact that there are no perfect metals. 
This fact can easily be reconciled with the (nearly) free 
electron model by the following arguments due to Drude 
and Sommerfeld. One assumes that the electrons in a 
metal are accelerated by an electric field for times smaller 
than the average time r needed to scatter off an impurity. 
As a result, the distribution of electron velocity reaches 
a steady state that is characterized by the finite Drude- 
Sommerfeld conductivity 



ne^T 
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where e, m, and n are the electron charge, mass, and 
density, respectively. This argument is essentially classi- 
cal. The assumption underlying Eq. (Q) is that electrons 
are particles, not waves. It took two pioneering papers, 
written by Dyson in 1953i and Anderson in 1958^ to call 
into question this assumption on theoretical grounds. 

Anderson's 1958 paper dealt with a quantum par- 
ticle hopping between neighboring sites on a three- 
dimensional lattice. Anderson added on-site potentials 
to each lattice site, which were taken to be random^Si with 
the standard deviation of the random on-site potentials 
far exceeding the hopping amplitude t. For this case, An- 
derson showed that any time-dependent wave packet that 
is initially localized both in space and in energy almost 
certainly remains so after an infinite time-evolution. An- 
derson further showed that the wavefunction envelopes 
decay exponentially fast away from their maxima, a phe- 
nomenon now called Anderson localization. Disorder was 



thus shown to play a role that far exceeds the one that it 
plays in the Drude-Sommerfeld picture. Disorder alone 
can turn the nearly free electron gas (the metallic state) 
into an insulator where wavefunctions are localized. This 
(Anderson) insulating state is completely different in na- 
ture from the band-gap insulator. 

In the late 1960s, Mott argued that, if the random- 
ness of the on-site potential is not too strong, a mobil- 
ity edge should separate delocalized states with energies 
close to the band center from localized states near the 
band edges Upon increasing the disorder strength, the 
two mobility edges approach each other until they meet 
at a critical value of the disorder strength. The possibil- 
ity raised by Mott that tuning the Fermi energy through 
the mobility edge triggers a transition from the metal- 
lic to an insulating state solely due to weak disorder in 
three dimensions raised the issue of the nature of such a 
transition. 

The decisive argument in support of interpreting the 
Anderson metal-insulator transition as a continuous tran- 
sition amenable to scaling ideas developed for critical 
phenomena in equilibrium phase transitions came from 
the scaling theory of Abrahams et alA Whereas localiza- 
tion occurs in three dimensions only for strong enough 
disorder, in one and two-dimensional versions of the An- 
derson model, eigenfunctions are always localized, irre- 
spective of the degree of the randomness of the on-site 
potentialsi^ In one dimension, the localization length ^ is 
as small as the mean free path L In two dimensions, ^ is 
typically much larger than i. 

Dyson's earlier 1953 paper dealt with a quantum par- 
ticle hopping along a one-dimensional lattice (if we use 
the language of the Anderson model). Dyson, however, 
considered the case that there are no on-site potentials. 
Instead, he took the hopping amplitudes ij.j+i = ij+i.j 
between neighboring sites j and j H- 1 to be random vari- 
ables. We refer to Dyson's model as the "random hopping 
chain" . The random variations of tj j-i-i cause scattering 
between Bloch waves, which is described by means of a 
mean free time r. Dyson calculated the density of states 
(DOS) per unit length p in the thermodynamic limit and 
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found that it diverges near the band center e = 0, 

= I ?i1 II (l^^l « (2) 
jerln |er| | 

Here, Pq is the DOS in the absence of randomness, which 
is finite and nearly constant near e = 0. In contrast, there 
is no divergence in the DOS, which is roughly equal to 
Pa{s), in Anderson's model with random on-site poten- 
tials. Further, it was found that the conductance g of the 
random hopping chain of length L is not exponentially 
small (as is the case in Anderson's model), but that it has 
an anomalously wide distribution precisely at e = 0, with 
average (g) ^ {£/L)-^/^. Again, the disorder has far more 
profoundly changed the properties of the chain than one 
would anticipate based on the Drude-Sommerfeld model. 

The origin of the difference between the random hop- 
ping chain at e = and Anderson's model with random 
on-site potentials is a symmetry that is present in the 
former, but absent in the latter. In the random hopping 
chain, for each eigcnfunction ipj with energy e, (—lyipj 
is also an eigenfunction, but with energy — e, and hence 
the energy eigenvalues always appear in pairs, ±e. This 
symmetry property follows from the fact that the Hamil- 
tonian of the random hopping chain only has matrix ele- 
ments that connect sites belonging to two distinct sublat- 
tices. Under this "sublattice symmetry" (SLS), the band 
center e = is special, and wavefunctions at that energy 
have anomalous localization properties. 

The band center of the random hopping chain can 
be viewed as an example of two mobility edges merged 
together i& Thus, it can serve as a relatively simple model 
to study the metal-insulator transition in the Anderson 
model in higher dimensions. The random hopping chain 
has a diverging localization length ^(e) ~ £\ In |er|| upon 
approaching the band center e = In the parlance of 
critical phenomena, it is convention to call the exponent z 
entering the scaling relation L ^ E^^^^ between a length 
scale L and an energy scale E the dynamical scaling ex- 
ponent, in which case the scaling relation ^(e) ^ £\ In lerj | 
corresponds to z = oo and signals an "infinite disorder 
fixed point" A 

Whenever the (typical) localization length ^ is much 
larger than the mean free path £, universality is expected. 
By universality is meant that transport and thermody- 
namic characteristics of a disordered sample only depend 
on dimensionality and some intrinsic properties such as 
the presence or absence of certain symmetries (an exam- 
ple is the SLS, but also time-reversal symmetry (TRS) 
and spin-rotation symmetry (SRS) play a role). Micro- 
scopic details such as the band parameters and the de- 
tailed form of the disorder should be irrelevant. For ex- 
ample, the Dyson singularity is universal in that it is in- 
dependent of the form of the probability distribution of 
the hopping randomness (as long as the standard devia- 
tion of the logarithm of the hopping amplitude is finite) . 

The outstanding open question in the problem of An- 
derson localization is to compute the evolution of the full 
probability distribution of physical observables as a func- 



tion of the disorder strength, and in particular to com- 
pute scaling exponents at the metal-insulator transition 
in two and three dimensions. New insights on disorder 
induced critical behavior in the problem of Anderson lo- 
calization has occurred on two fronts that we shall review 
below. The first front stems from the renewed interest 
in the random hopping problem in samples of dimen- 
sions larger than one. and in its three incarnations asso- 
ciated to TRS and SRS. The second front stems from the 
introduction by Altland and Zirnbauer of four symme- 
try classes associated to dirty superconductorsiS In both 
cases, the existence of disorder induced critical points was 
established in a quasi-one dimensional geometry (sample 
with a length L much larger than its width) and shown 
to be identical to the Dyson critical point of the ran- 
dom hopping chainpi2ili Progress also took place in two 
dimensions with the exact computation of some critical 
exponents among these additional symmetry classes of 
Anderson localizationiiS We shall review these develop- 
ments in the following sections. 



II. TRANSFER MATRIX ANALYSIS IN 
QUASI-ONE DIMENSION 

Most detailed theoretical results for the localization 
problem have been obtained in one-dimensional and quasi 
one-dimensional disordered systems. The width of a 
quasi one-dimensional conductor is defined by the num- 
ber of propagating channels N at the Fermi level. A 
one-dimensional conductor has = 1. 

The localization problem in (quasi) one dimension 
has been tackled using a variety of theoretical methods, 
most notably field theory, network models, and trans- 
fer matrix approaches. In the field-theoretical approach 
the localization problem is mapped onto a field the- 
ory known as the (one-dimensional) "non-linear sigma 
model" (NLSM). The NLSM has been the canonical 
framework to establish universality, as the same field 
theory is obtained for a wide class of microscopic dis- 
order. An essential step in the derivation of one univer- 
sal field theory is a separation of length scales (in this 
case i -€1 ^ = 0{N£)) which restricts the applicability of 
this method to quasi one-dimensional systems with many 
channels, ^ 1. The other two methods do not require 
the limit A^ ^ 1. However, equivalence with the NLSM 
and. thus, universality, is to be expected for A^ ^ 1 only. 

All three approaches have contributed to the present 
understanding of the localization problem. The field- 
theoretical approach suggested the existence of critical 
points induced by disorder belonging to different uni- 
versality classesii^ It has also played an important role 
in the characterization of the mesoscopic fiuctuations^ 
The mapping onto a network model allowed the numer- 
ical computation of the correlation (localization) length 
scaling exponent in the plateau transition corresponding 
to the lowest Landau level in the integer quantum Hall ef- 
fect (IQHE)^ and has also been applied to the new uni- 
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versality classesii^ The transfer matrix approach was im- 
plemented numericaUy to compute the scahng exponents 
of the Anderson locahzation in three dimensionsii2*i£ It 
also yielded important analytical results for the quasi 
one-dimensional case, which we now review. 

Let us first illustrate how the transfer matrix approach 
is applied to the random hopping problem in one dimen- 
sion {N = 1). We ignore the electron spin and assume 
that disorder is weak compared with the band width. 
Since we are interested in properties close to the band 
center £ = 0, the energy spectrum {e{k) oc — cos/c) can 
be linearized about it (i.e., about k ~ ±7r/2) and the 
wavefunction ipj can be written as 

i;j^e'^^My) + e-'^'My) (3) 

with y = jx (lattice spacing) , and ip,. and ip-^ correspond- 
ing to waves moving to the right and to the left, respec- 
tively, at velocity Vp. The continuum limit of the Hamil- 
tonian for the random hopping chain can then be rep- 
resented by the stationary Schrodinger equation (h = 1) 
for a spinor = [ipj., ip-y )"^, 

(H-£)* = 0, n^-ivpT^dy-Viy), (4) 

where 2 3 S'^'^ the Pauli matrices in the left- 
mover/right-mover grading. From Eq. (O we infer that 
multiplying ipj by the factor (—I)-' = e"^^ amounts to 
exchanging ip^ and ipx, or equivalently, to the transforma- 
tion \1/ Ti^ . Hence the SLS of the Hamiltonian, that 
an eigenfunction ipj with energy e has its partner {—lyipj 
with energy — e, reads, in the continuum language, 

tiHti ^ ~n, (5) 

which implies that the scattering potential V can have 
terms proportional to T2 and T3 only. Without SLS, V 
also has terms proportional to ti and the 2x2 unit ma- 
trix, which do not satisfy the symmetry jSJl- 

The transfer matrix A4^{L,0) relates the wavcfunc- 
tions \E'(0) and ^'(-L), 

*(L) = A1^(L,0)*(0). (6) 

Knowledge of A4 is sufficient to calculate the Landauer 
conductance g: the eigenvalues of the product M^Ml 
come in the pair exp(±2a;), which is related to g by 

Transfer matrices satisfy a multiplicative rule: the trans- 
fer matrix of a disordered segment of length Li + L2 is 
simply the product of the transfer matrices of the disor- 
dered segments of lengths Li and L2; 

Me{Li + L2, 0) = M,{Li + L2, Li)Me{Li,0). (8) 

Using the Born approximation to calculate the transfer 
matrix of a segment of length much smaller than the 
mean free path £, the transfer matrix of the full system 




FIG. 1: The "radial coordinate" of the transfer matrix makes 
a Brownian motion on an associated symmetric space. 



is then found by repeated matrix multiplications. What 
makes an analytical solution possible is that each multi- 
plication changes the transfer matrix only slightly. As we 
increase the system size from L to L + SL, the random 
potential V{y) in the interval y G [L,L + SL] changes 
the transfer matrix from Mg{L, 0) into Me{L + SL, 0) = 
Me{L + 6L, L)Mg{L,0). Drawing an analogy with the 
Brownian motion, we may regard L as a time and the 
random potential V{y) as a random force; A4^{L, 0) then 
performs, as a function of length L, a "random walk" in 
the manifold of allowed transfer matrices. 

The structure of the symmetric space is determined by 
the symmetries of the transfer matrix A4^. These follow 
from the fundamental symmetries of the Hamiltonian Ti.: 
Hermiticity, TRS, and the SLS. For the transfer matrix, 
these imply 

MI{L)t^M,{L)^t^, (9a) 
TiM:{L)ti=M,{L), (9b) 
T,M+,iL)T,=M_,iL), (9c) 

respectively. Here we have used the shorthand notation 
Ale(L) := A4e{L, 0). These conditions define a manifold 
(Lie group) on which the transfer matrices A4e live. 

In view of Eq. Q , we are interested in the radial diffu- 
sion only, related to the eigenvalues exp(±2a;) of A^^A^J. 
Mathematically, this random walk is reformulated into 
a random walk on an associated "symmetric space" 
which is obtained after identifying transfer matrices 
that have the same product M^Ml (Fig.^. Symmetric 
spaces have a well-defined natural metric, and if the ran- 
dom potential has a Gaussian distribution without long- 
range correlations, one finds that the random walk can 
be described by the diffusion equation for the symmetric 
space ifi* 

The SLS (|9c|l poses an extra condition (commutativity 
with Ti) on the transfer matrix at £ = 0, that is absent for 
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transfer matrices at finite energy or for transfer matrices 
of a wire without SLS. In other words, e breaks the extra 
symmetry of A4e=o, which we also call SLS. This fun- 
damentally changes the structure of the symmetric space 
and, hence, the solution of the localization problem: at fi- 
nite energy e, the radial coordinate x {i.e., the logarithm 
of the eigenvalue of M^Ml) performs a biased random 
walk, resulting in an exponentially small conductance for 
large L: generically, all wavefunctions are localized in one 
dimension. On the other hand, at e = one finds that 
as X performs an unbiased random walk, the root-mean- 
square of the radial coordinate x grows proportionally 
to Vl whereas t he ty pical values of the conductance de- 
cays as exp(— a-\/L/^) where a is a positive constanti^i*^ 
Since there is a finite probability density cx L~^/^ to find 
X at the origin, the ensemble averaged conductance de- 
cays only algebraically, (g) cx L~^^^^ 

A quasi one-dimensional system with spinful electrons 
is obtained by replacing the wavefunctions t/'j. and ipi by 
2A^-component vectors. The disorder potential V and the 
transfer matrix Ai then become 4A^ x 4A^ matrices. The 
eigenvalues of M come in Z?-degenerate pairs exp(±2a;„), 
n — 1,2, . . . , N* , where the degeneracy D depends on 
the symmetries and N* ~ 2N/D. The Landauer con- 
ductance reads24 

AT* 

ff = ^E^;^- (10) 

; cosh Xn 

Upon increasing L, M. performs a random walk in the 
manifold of allowed transfer matrices. If the matrix ele- 
ments of the disorder potential have identical and inde- 
pendent Gaussian distribution without long-range cor- 
relations, this random walk again is described as dif- 
fusion on an associated symmetric space. The corre- 
sponding diffusion equation (or Fokker-Planck equation) 
for the probability density of the radial coordinates Xn, 
n = 1, . . . , N* reads 

ap _ j_ ^ _a_ 5 1 

where J is a Jacobian describing the transformation to 
radial coordinates. It has the functional form 

J = Y[ (2a;,) Yl Yl sinh""* {x, ± Xk), (12) 

j k<j ± 

where m; and mo± are numerical constants describing the 
geometric structure of the symmetric space (see TablePl. 
Finally, 7 = (mo+ + mo-){N* — l)/2-|- 1 -|-to/ is another 
numerical constant. 

To determine the appropriate symmetric space, one be- 
gins from the three symmetries listed above (TRS, SRS, 
and SLS). A statistical ensemble of Hamiltonian is clas- 
sified into symmetry classes according to its invariance 
properties under the transformations induced by these 
symmetries. This gives 6 symmetry classes with the cor- 
responding symmetric spaces given in the first 6 rows 



from Table ^ In the context of quasi one-dimensional 
localization for the standard symmetry classes (O, U, S 
m Table D), Eq. ^ was first derived by Dorokhov^ 
It was re-derived independently by Mello, Pereyra, and 
Kumar)^ and is now called the DMPK equation^ The 
three symmetry classes with SLS are called chiral classes. 
Of the three chiral random matrix theories that describe 
the energy level statistics on the scale of the mean- 
level spacing of the random hopping problem^S^ the chU 
class has also been applied to the chiral phase transition 
in quantum chromodynamics^ The diffusion equation 
associated to for the localization problem with 

SLS was derived by Altland, Simons, and two of the 
authors 

Four more symmetry classes are obtainedSi when con- 
sidering the effect of disorder on Bogoliubov quasipar- 
ticles in a superconductor described at the mean-field 
level^ In a superconductor, the quasiparticlc wavefunc- 
tions at excitation energy e are cigenfunctions of the 
Bogoliubov-de Gennes (BdG) equation, 

(7^-e)* = 0, W= (^_^, (13) 

which has an additional 2x2 grading, the particle- 
hole grading, when compared to the Schrodingcr equa- 
tion. Here h (A) is a 4A^ x \N Hermitian (antisym- 
metric) matrix representing a single-particle Hamiltonian 
(superconducting order parameter). The Hamiltonian 
H. in the BdG equation satisfies the symmetry relation 
7i7i*7i = —Ji (particle-hole symmetry), where 71 is the 
Pauli matrix in the particle-hole grading. For the trans- 
fer matrix particle-hole symmetry implies 

7iX+,(i)7i=Al-.(^). (14) 

At e = 0, the symmetry requirement (|14|) modifies the 
structure of the manifold of allowed transfer matrices and 
thus qualitatively modifies the localization behaviori^i 
The presence of particle-hole symmetry gives four more 
symmetry classes, referred to as C, CI, D, and Dili, de- 
pending on the presence or absence of TRS and SRS, 
see Table ID At the level of mean-field theory, charge is 
not a good quantum number. Only energy is conserved 
as is spin in the presence of SRS. Anderson localization 
applied to BdG quasiparticles thus aims to describe the 
statistics of the global and local DOS of BdG quasiparti- 
cles and the thermal or spin conductance. The diffusion 
equation (|11() for the four BdG symmetry classes was re- 
ported at the end of the 20th centurypii 

In total, there are 10 different symmetry classes, yield- 
ing 10 symmetric spaces for the transfer matrix ensemble. 
For four of those symmetry classes (U, chU, CI, Dili), 
an exact solution of the diffusion equation Hll|l is possi- 
ble. For the remaining classes, asymptotic solutions for 
L <C T^l and L ^ are known. We first review the 
asymptotic for small and large lengths, and then discuss 
the exact solutions. 

In the regime i <^ L <^ Nl where transport is dif- 
fusive, moments of the conductance g can be obtained 
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TABLE I: Classification of symmetry classes for disordered quantum wires. Symmetry classes are defined by the presence or 
absence of time-reversal symmetry (TRS) and spin-rotation symmetry (SRS), and by the further symmetries of the disorder: 
sublattice symmetry (SLS) (random hopping model at the band center), and particle- hole symmetry (PHS) (zero-energy 
quasiparticles in superconductors). For historical reasons, the first three rows of the table are referred to as the orthogonal (O), 
unitary (U), and symplectic (S) symmetry classes when the disorder is generic. The prefix "ch" that stands for chiral is added 
when the disorder respects a SLS as in the next three rows. Finally, the last four rows correspond to dirty superconductors 
and are named after the symmetric spaces associated to their Hamiltonians. The table lists the multiplicities of the ordinary 
and long roots mo± and m; of the symmetric spaces associated with the transfer matrix. Except for the three chiral classes, 
one has mo+ = nio- = mo- For the chiral classes, one has mo+ — 0, mo- = mo. The table also lists the degeneracy D of 
the transfer matrix eigenvalues, as well as the symbols for the symmetric spaces associated to the transfer matrix M and the 
Hamiltonian Ti. The last three columns list theoretical results for the weak-localization correction Sg for £ <^ L <ti N£, the 
average of Ing at L ^ N£, and the DOS near e = 0. The results for {Ing) and p(e) in the chiral classes refer to the case of A'' 
even. For odd N, {\ng) and p{e) are the same as in class D. 
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as a power series in L/Nim^^^ii^ The leading term in 
the expansion for the average (17) is the Drude conduc- 
tance (cx N£/L), cf. Eq. Q). The leading correction to 
the Drude conductance is known as the weak locahzation 
correction which, unUke the Drude conductance, depends 
on symmetry (see Table Pi. When negative, 6g is consid- 
ered as a precursor to exponential localization. For the 
unitary and the three chiral symmetry classes, Sg = 0. In 
the BdG symmetry classes, the sign of Sg depends on the 
absence or presence of SRS, a situation similar to the one 
for the orthogonal and symplectic cases in the standard 
symmetry classes, but the behavior at L ^ N£ is very 
different from the standard classes. 

The probability distribution of the conductance be- 
comes very broad if L 3> N£ and is poorly characterized 
by its mean. For the standard symmetry classes, the 
chiral classes with an even number of channels, and for 
the BdG symmetric classes CI and C, In 5 becomes self 
averaging and thus represents the distribution well. In 
these classes. In 17 increases linearly with L, correspond- 
ing to exponential localization. We refer to Table ^ for 
the relevant results. For the remaining cases, In g scales 
proportional to (L/iV£)^/^ (implying the diverging local- 
ization length), and its fluctuations are of the same order 
as its mean. This is a signature of anomalous localiza- 
tion that reflects itself with the anomalously slow decay 
~ y/WI/L of the mean conductance. 

The exact solution of the diffusion equation Ijllfl is pos- 
sible if the multiplicity of the ordinary root nio is two. 
The exact solution has been used to find the first and sec- 
ond moments of the conductance for all lengths. For the 
chU, CI, and DIII symmetry classes, the first moment is 



given bjiii*^2i2i 

00 

(■9)- = i- + i^(-l)"(^+i)e--="V8.^ (15) 

71=1 

s 3 ■^-^ \ s TT^n^ / 

n— 1 ^ ^ 

19°° 1 

idUu = - + ^-4Ee"^'"'/'^4^' (17) 
s 3 TT^n^ 

n— 1 

respectively, where s = L/AN£, in the universal limit 
N 00. The exact result shows explicitly that localiza- 
tion is non-perturbative in L, explaining why only very 
little progress can be made with standard perturbative 
techniques. 

In Fig. 121 we show results of numerical simulations of 
the mean and variance of the conductance for the so- 
called quasi-lD random flux model where tij = e'^'-^ 
with random 9i j and which has the SLS but no TRS, 
together with the exact results with ^ as a fitting pa- 
rameter. What is specific to the chiral symmetry classes 
realized at e = is an even-odd effect in N: as in the ran- 
dom hopping chain, (g) decreases algebraically as L~^^^ 
with L for odd N whereas it decays exponentially fast 
with L for even A^. This remarkable even-odd effect is 
reminiscent of the one for the spin gap in the number of 
legs of S" = 1/2 spin ladder systemsi^ For the diffusive 
regime (^ ^ L <C N£), the second cumulant of g is called 
the universal conductance fluctuations, and its value in 
the chiral classes is twice as large as the one in the cor- 
responding standard symmetry classes. The numerical 
results at £ ^ show the behavior of the unitary class as 
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FIG. 2: The system size (L) dependence of (a) mean {{g}) and 
(b) variance {{g'^) — {g)^) of the conductance in the random 
flux model. The numerical results at the band center e — 
are denoted by o (iV = 15) and □ (iV = 16), and the 
exact solution (at ^ 1) from the diffusion equation are 
shown by full (A*': odd) and dotted (A'^: even) lines. For the 
variance, the result off the band center e = 0.1 with = 32 
is also displayed by A for comparison. It is well-fitted by the 
analytical result (broken line) of the unitary symmetric class 
which is obtained by solving the NLSM exactly in quasi-one 
dimensions^ 



N£, the DOS is constant and given by pg = ND/{-KVp) 
for all symmetry classes. This is also the behavior of the 
DOS for the standard symmetry classes when \eTc\ < 1. 
For small energies, \stc\ < 1, the remaining symmetry 
classes, for which the energy £ = plays a special role, 
arc much more interesting. The expectation that level 
repulsion would enforce a vanishing of the DOS captured 
by random matrix theory is wrong whenever the long root 
m; vanishes. As shown in tabled p(e) shows a power law 
behavior with logarithmic corrections whose exponent for 
the chiral classes with N even is not captured by the 
chiral random matrix theory. For the chiral classes with 
N odd and the BdG classes D and Dili, the DOS exhibits 
the Dyson singularity. 

The diffusion equation (|ll|l is a beautiful example of 
one-parameter scaling. Here we'd like to exert one word 
of caution about the role played by symmetry: our for- 
mal discussion was conditional on the assumption that 
all matrix elements of the weak potential V have iden- 
tical, independent, and Gaussian distributions. If these 
conditions are not met, the conductance and DOS dis- 
tributions will depend on microscopic details. Obviously, 
in such a case, symmetry need not be the only factor 
in determining the conductance and DOS distributions. 
In particular, if symmetries are only weakly broken (for 
example, when TRS is broken by a weak magnetic field 
or when a small e > breaks the SLS of Al^), one 
may end up in a cross-over between two or more of the 
symmetry classes listed in Tabled However, in the limit 
of large channel number A^, universality is believed to 
be restored. Symmetry is the only player in the limit 
N ^ 1 (keeping L/N£ fixed, to ensure a fixed position 
in the crossover to localization), so that the notion of a 
symmetry class can be elevated to that of a universality 
class. 



a result of the breaking of the SLS 

What is the relation between the divergence of the 
localization length and the DOS singularity originally 
discovered by Dyson? In order to see this, one has to 
adapt the transfer matrix approach to the calculation of 
scattering phase shifts. Phase shifts determine the DOS 
through the Friedel sum rule. Interestingly, information 
about the scattering phase shifts can be obtained from 
the standard transfer matrix approach if we reconsider 
the transport problem in the presence of an imaginary en- 
ergy e = iujm^ Unlike a real energy, an imaginary energy 
does not change the symmetry of the transport problem. 
The energy iuj simply yields a drift term in the diffu- 
sion equation Hll|) . but otherwise allows the transport 
problem to be solved by essentially the same methods 
as outlined previously. Hence, just like the conductance, 
the DOS is thus completely determined by the geometry 
of the corresponding symmetric spaces! 

For energies of the order of the inverse of the time 
Tc ^ N'^£/vp needed to diffuse through a wire of length 



III. PROGRESS IN TWO DIMENSIONS 

One of the best studied disorder-induced critical points 
in two dimensions is the plateau transition in the IQHE. 
Upon changing the filling fraction of the (lowest) Landau 
levels by applying a magnetic field, say, the localization 
length for states at the Fermi level diverges at some crit- 
ical value B^, 

^^\B-B,r, (18) 

at which the Hall conductance jumps by one in the unit 
of e^/h. Here, the exponent for the localization length 
ly is known numerically to be 2.3 < v < 2.4 within the 
framework of Anderson localization^ The single-particle 
wavcfunction 'ip(T) the plateau transition is multifrac- 
tal, i.e., the so-called inverse participation ratio Vq scales 
with the linear size L of the system as a power law with 
the non-linear function of scaling exponents T{q), 

Vg^ f d2r|^(r)|"« (19) 
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The spectrum of exponents T{q) is beheved to be univer- 
sal. 

The successful theory of the plateau transition in the 
IQHE must predict u and T{q). Decrypting the critical 
(conformal field?) theory describing the plateau transi- 
tion remains a tantalizing open problem. On the other 
hand, there has been some progress on two-dimensional 
criticality in the chiral and BdG universality classes. In 
this section, we will review these results with main em- 
phasis on exact results. 

Within the two-dimensional chiral symmetry class, the 
exact computation of the multifractal scaling exponents 
r(g) was carried out for the Hamiltonian describing a par- 
ticle with relativistic dispersion in a random white-noise 
correlated vector potential of vanishing mean A^{r), 

fj.= l,2 

(a^ is a Pauli matrix). This problem of Anderson local- 
ization is a fine-tuned model with symmetry of the chU 
type that realizes a line of critical points as a function 
of the variance gA of A^i^ The band center e = is a 
mobility edge, where a critical state is located, as was 
the case in the random hopping chain; see Fig. 3. The 
exact computation of the T((jf)-spectrum for this critical 
wavefunction is possible by exploiting similarities to the 
problem of directed polymers in random medial For 
|g| < 27T /gA , T{q) is a quadratic function of q whereas 
it is linear in q outside this region (Fig. 3). A similar 
behavior was also observed numerically for critical wave- 
functions in the IQHEi^ The DOS too can be computed 
exactly close to the band center. Its dependence on en- 
ergy e is algebraic, p(e) ^ \e\^. The exponent /3 also dis- 
plays a non-analytic dependence on g^mi^^ Both singu- 
larities can be ascribed to the non-analyticity at 5^ — 27r 
of the dynamical exponent^Siii 

(1 + ^, for5^<27r, 
z{gA) = I (21) 
IVit-l' for5A>27r. 

This non-analyticity as a function of can be inter- 
preted as a freezing transition by analogy to the non- 
analytic dependence on temperature of the free energy 
in the random-energy model for spin glasses4ii 

The random vector potential problem H2()|l is the sim- 
plest two-dimensional continuum model with randomness 
that encodes the chiral symmetry. For more generic mod- 
els with SLS, the critical line is unstable to marginal 
perturbations (in the renormalization group sense) com- 
patible with the chiral symmetry, which renormalize 5^ 
to strong coupling..^ Gade introduced a NLSM that de- 
scribes generic two-dimensional random hopping prob- 
lems with SLS, for which she found the diverging DOS 

p{e) 1- expire] \ne\'^) (22) 




FIG. 3: (a) Typical wavefunction amplitude |i/;(r)p and (b) 
its exact r(q)-spectruni (solid line) for the random Hamilto- 
nian i'AOl at — 1.2. 



sufficiently close to the band center. Here, c is a non- 
universal constant and k = 1/2 according to Gadeii^ 
This result can be derived by dimensional analysis if 
one assumes that energy and length scales are related by 
Ine ~ — z(L)lnL with z{L) ~ InL. However, Motrunich 
et al. conjectured that k = 1/2 is only a transient as 
the true asymptotic dependence of the effective dynam- 
ical exponent on L is z{L) ^ \/hiL, as would follow if 
the L dependence of z{L) was given by Eq. H21|) with 
an effective running coupling constant ^ InL.*-^ If so, 
the scaling exponent k = 2/3 follows instead of k = 1/2. 
The latter scaling analysis was confirmed by a functional 
renormalization group calculation for a two-dimensional 
random hopping problem in Ref . l43i 

According to the theorem of Anderson on the effect of 
disorder in an s-wave BCS superconductor, non-magnetic 
impurities are largely an irrelevant effect. In an uncon- 
ventional BCS superconductor, here defined by an or- 
der parameter with a nonzero angular momentum, non- 
magnetic impurities have a much stronger effect. With 
the realizations of unconventional superconductivity such 
as high-Tc superconductors, spin triplet superconduc- 
tivity in ruthenates, two-dimensional organic materials, 
paired states in the fractional quantum Hall effects etc., 
searching for non-pcrturbative effects in dirty, unconven- 
tional superconductors gained in urgency in the mid 90s. 
With SRS and no strong breaking of TRS, BdG quasi- 
particles are localized in two dimensions. In this case the 
DOS is depressed close to the band center as expected 
from level repulsion, p{e) ^ e with TRS (CI class) and 
p{e) ^ without TRS (C class) On the other hand, 
as is the case with the IQHE, strong breaking of TRS 
opens up the possibility of a new critical point induced 
by disorder, the thermal (spin) Hall transitions^ Using 
a network model representation of the spin Hall tran- 
sition in Class C Gruzberg et al. reduced the problem 
of computing the exponent v for the diverging localiza- 
tion length and the exponent f3 for the power law decay 
of the DOS to the classical problem of bond percola- 
tion. They found v = 4/3 and (3 = 1/7 at the plateau 
transitionpi^ No results of that magnitude are presently 
available when SRS is broken in a dirty BdG but the exis- 
tence of a plateau transition for the Hall thermal current 
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has been established numericaUy together with a very 
rich phase diagram 4^ 

IV. SUMMARY AND OPEN PROBLEMS 

Recent theoretical progress has established that criti- 
cal phenomena at mobility edges in quasi-one dimension 
arc of the Dyson type. Their properties arc much bet- 
ter understood than with Anderson localization in two 
dimensions in all but few cases. In this article, we have 
reviewed two-dimensional models for which the multi- 
fractal spectrum or some critical indices can be obtained 
exactly. Within a classification based on symmetries of 
disorder induced critical points in the problem of Ander- 
son localization, the chiral symmetry class appears to be 
the simplest both in quasi one and two dimensions. 

The physics of the chiral symmetry classes is closely 
related to the physics of the pseudogap in Pcicrls sys- 
tems and to the physics of the chiral phase transition 



in quantum chromodynamics, which we could not dis- 
cuss in this limited space. Chiral and BdG classes are 
also related to classical and quantum random spin sys- 
tems. For examples, the Dyson singularity in the ran- 
dom hopping chain is connected to the random singlet 
phase in one-dimensional quantum spin systemsj^^ the 
two-dimensional chiral symmetry classes are related to 
the two-dimensional gauge glass modelj^^i^ and the clas- 
sical random bond Ising model in two dimensions can 
be mapped onto a two-dimensional network model be- 
longing to the symmetry class D. In the last example, 
this connection was used to demonstrate the multifrac- 
tal behavior— of the correlation function of a dual order 
parameter and a freezing phenomenoniSi 

Another important open issue is that of the competi- 
tion between disorder and interactions in the chiral and 
BdG universality classes. The situation seems to be un- 
der better control for the BdG symmetry classes than in 
the standard symmetry classes— 
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